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BACKGROUND OF THE INVENTION 
1* Field of the Invention 

The present invention relates generally to registration of images and, 
more particularly, to a method for blind cross-spectral image registration. 

2. Prior Art 

Image registration is the process of aligning two images of the same 
scene so that corresponding points in the scene are placed in identical pixel positions. 
Standard full-color reproductions use precisely registered images for each of the 
component colors. Similarly, false color images combine registered image planes 
from various spectra to reveal important details not readily apparent in the individual 
images. For remote sensing, registration of infrared to visible spectra is especially 
important for measuring vegetation, detecting ocean currents, and tracking hot spots in 
forest fires. Registration of images taken at different times is typically used to 
identify changes between the images. 

The prior art for the problem of image registration generally falls into 
two different approaches— feature-based and blind. Feature-based registration 
attempts to identify edges, comer points, contours, or other features that are common 
to two images, and then uses standard geometric transforms to compute the mapping 
between the pairs. The problem of identifying those features is complicated by the 
fact that edge features in infrared images are related to temperature variations, and 
these edges may not be present in the visible spectrum. Likewise, some features in the 
visible spectrum may disappear in the infrared spectrum. Consequently, feature-based 
registration is mainly concerned with locating features common to both images, and 
rejecting features that are exclusive to one image. The problem becomes difficult 



when relatively few features are common between the images. For example, a pair of 
aerial images of an agricultural region may show relatively uniform intensity in the 
visible spectrum, and highly textured intensity in the infrared spectrum. Each feature 
evident in the visible spectrum may map to many possible candidates in the infrared 
5 image. 

The second approach to the problem is to register images blindly by 
maximizing some criterion that depends on the quality of a candidate registration. The 
second approach completely avoids the problem of finding a subset of features 
common to both images, and matching the features to each other. Typical criteria for 

1 0 bUnd registration are to minimize the sum of squared differences of pixel values or to 
maximize the normalized correlations of the images. Perhaps the most powerful 
criterion is the maximization of mutual information which is particularly effective 
when one image differs from the other in a rather complex way, such as might be 
observed due to changes in the illumination source position, image modality (X-ray 

1 5 and MRI), or spectral channel (visible and infrared). It has been used effectively in 
practice to register PET, MR, and CT medical images. 

A major potential disadvantage of mutual-information-based methods 
is the large computational overhead required to compute the joint distributions 
between two images for many different relative alignments of the images. To 
20 overcome this disadvantage, those in the art describe nonlmear iterative methods that 
reduces substantially the number of different relative aUgnments that need to be 
exammed. Although the non-linear iterative methods use a sum of square differences 
of pixel values as the criterion for registration quality, it is known to use mutual- 
information criterion in its place. 

25 Although the non-linear iterative methods, like all blind-regisfration 

algorithms, avoid the cost of identifying correspondmg features, the computation is 
expensive, even in the iterative form of the method. For each relative position of the 
images considered, a joint distribution of pixel values needs to be computed, which 
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involves a number of operations proportional to the size of the image. Coarse-to-fine 
techniques known in the art help reduce this cost. Nevertheless, the algorithm must 
examine several different displacements at maximum detail and many more at lesser 
detail, and each examination involves access to all of the pixel values at that level of 
detail. 

SUMMARY OF THE INVENTION 

Therefore it is an object of the present invention to provide a method 
for registration of images with comparable quality as methods which employ 
maximization of mutual information but with lower computational complexity. 

The registration methods of the present invention concern a fast 
technique for registering image pairs from visible and infrared spectra that differ by 
translation, small rotations, and small changes of scale. The main result of the 
registration methods of the present invention is a nonlinear prefiltering and 
thresholding technique that substantially enhances the cross-spectral correlation, 
provided that the image pairs have many features in common. The non-linear 
prefiltering and thresholding techniques provided are used in conjunction with a 
Fourier-based normalized correlation method to perform fast cross-spectral 
registrations. In the absence of such prefiltering, local reversals of contrast from 
image to image tend to impair the quality of correlation-based registrations. 

The registration methods of the present invention are blind in that they 
do not identify specific features in both images to use for alignment. Instead, they 
compute the translation that maximizes the overall normalized correlation of the 
filtered images. Small rotations and scale changes can be recovered by computing the 
translation displacement in several different regions of the image pairs. Fourier 
techniques for computing normalized correlations greatly reduce computational costs, 
and eliminate the necessity to use iterative search techniques to hold computation 
costs down. 



Accordingly, a method for registration of first and second images out of 
registration is provided. The method comprises the steps of: (a) making the edges in 
the first and second images more prominent; (b) thresholding the first and second 
images fi-om the previous step using a threshold for which N percent of the pixels of 
each of the first and second images are over the threshold; (c) reducing the resolution 
of the first and second images fi-om the previous step; and (d) registering the first and 
second images of reduced resolution fi-om the previous step. 

Preferably, the method further comprises the step of blurring the first 
and second images from the thresholding step. The blurring step preferably comprises 
filtering each of the first and second images firom the thresholding step such that each 
pixel therein is thickened by a predetermined number of pixels in a square array that 
extends the predetermined number of pixels in all four directions from a central pixel 
The method preferably also fiirther comprises the step of increasing the resolution of 
the registered first and second images firom the registering step. 

Also provided are a computer program product for carrying out the 
methods of the present invention and a program storage device for the storage of the 
computer program product therein. 

BRIEF DESCRIPTION OF THE DRAWINGS 

These and other features, aspects, and advantages of the methods of the 
present invention will become better understood with regard to the following 
description, appended claims, and accompanying drawings where: 

Figures la and lb illustrate a wetiands image, with Figure la 
illustrating the red channel intensity of the wetiands image and Figure lb illustrating 
the infirared channel intensity. 

Figures 2a and 2b illustrate the images of Figures la and lb, 
respectively, after filtering to enhance edges using a filter coefficient of C- 8.5. 



Figures 3a and 3b illustrate the images of Figures 2a and 2b, 
respectively, after thresholding to create 80% white pixels. 

Figures 4a and 4b illustrate the images of Figures 3a and 3b, 
respectively, after thickening by 5. 

5 Figures 5a and 5b illustrate the images of Figures 3a and 3b, 

respectively, after thickening by 9. 

Figures 6a and 6b illustrate the images of Figures 3a and 3b, 
respectively, after thickening by 17. 

Figure 7a illustrates the image of Figure 3a after resolution reduction 

10 by 64. 

Figure 7b illustrates the image of Figure 6a after resolution reduction 

by 64. 

Figures 8a and 8b illustrate an agricultural image, with Figure 8a 
illustrating the red channel intensity of the agricultural image and Figure 8b 
1 5 illustrating the infi-ared channel intensity. 

Figures 9a and 9b illustrate a forestry image, with Figure 9a illustrating 
the red channel intensity of the forestry image and Figure 9b illustrating the infi*ared 
channel intensity. 

Figures 10a and 10b illustrate an urban image, with Figure 10a 
20 illustrating the red channel intensity of the urban image and Figure 1 Ob illustrating the 
infrared channel intensity. 

Figure 1 1 illustrates a graphical summary of Receiver-Operating 
Characteristics (ROC) data. 
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Figure 12a illustrates ROC curves for registrations of images captured 
in the red spectrum and images captured in the blue and green spectra. 

Figure 12b illustrates ROC curves for registrations between images 
captured in the red spectrum and images captured in the infrared spectrum. 

Figure 13 illustrates sensitivity of ROC curves to center coefficient of a 
9-point filter, where all filters except C ^ 8.0 are thresholded at 80% and the filter for 
C = 8.0 is thresholded at 50%. 

Figure 14 illustrates the sensitivity of ROC curves to the edge 

threshold. 

Figure 15 illustrates the sensitivity of ROC curves to edge-thickening 

coefficient. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

The registration methods of the present invention use an alternative 
strategy for cross-spectral image registration that takes advantage of Fourier 
techniques to reduce the complexity of normalized correlation. One can compute 
normalized correlations of two images for all relative integral displacements of the 
two images for a small constant times the cost of the normalized correlation at a single 
position. Since the cost for one normalized correlation is about the same cost as one 
mutual information evaluation at the same level of detail, it is clear that the Fourier- 
based normalized correlation method enjoys a computational advantage over the 
mutual-information methods of the prior art. The pixel-based normalized correlation 
registration methods whose timing is known in the art is much slower than its Fourier- 
based counterpart. 

The main difficulty in using fast correlation is that cross-spectral 
images typically have poor correlations. An objective of the present invention is to 
find a way to process images so that they can be registered accurately by means of fast 



normalized correlation. Those in the art have discussed the general problem of 
registering images across spectra and commented that global measures of registration 
accuracy usually work poorly in such cases. Any features that are exclusive to one 
image or another cause problems for normalized correlation, mutual-information, and 
5 other global measures because the exclusive features degrade the similarity measures. 
Those in the art processed the images to enhance image similarity, and then used local 
correlation rather than global correlation. They used an iterative scheme similar to 
that of the non-linear iterative methods to find a registration that maximizes the sum of 
the local correlations. The registration methods of the present invention also process 
10 the images, but does so in a way that enables fast global correlation to succeed. 

The main results of the registration methods of the present invention lie 
in the combination of an image preprocessing method and fast normalized correlation 
to register cross-spectral images with about the same quality as maximization of 
mutual information but with lower computational complexity. The preprocessing uses 
15 both edge enhancement and thresholding, with optional blurring (alternatively referred 
to in the art as thickening), which may be useful in conjunction with coarse-to-fine 
registration. For a moderate-sized data base, the registration methods of the present 
invention registered cross- spectral images about as well as mutual-information 
registration at full resolution, and was slightly inferior in quality at lower resolution. 

20 The registration methods of the present invention are based on 

normalized correlation of nonlinearly filtered images. Before expanding on the same, 
a general overview is first given. 

Figures la and lb illustrate the nature of the registration problem. The 
two images in Figures la and lb are aerial photos that are slightly out of registration. 
25 The image of Figure 1 a illustrates an image in the red visible spectrum, while the 
image of Figure lb is in the infrared spectrum. Note that there are some intensity 
inversions from image to image, but in many regions there is no intensity inversion. 
Registration techniques based on normalized correlation tend to perform poorly in 
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these circumstances. Normalized correlation measures how well an affine mapping of 
image intensities explains the differences between the intensities of two images. Local 
intensity inversions tend to fall outside an affine mapping, thereby lowering 
correlation values. 

5 When cross-spectral region intensities correlate poorly, it is reasonable 

to use region boundaries rather than intensities for registration purposes. Because of 
changes of intensity in different spectra, the boundaries visible in one image generally 
do not correspond to boundaries in the other image in a one-to-one fashion. 
Moreover, even where they correspond, the detected boundaries in the two images 
1 0 may have different pixel structure, and therefore may not register well This tends to 
reduce the correlation coefficient, and makes precise registration very difficult. 

A first step in the registration methods of the present invention is to 
filter the images with an edge- enhancement filter to make the edges prominent. 
Prefiltering is illustrated in Figures 2a, 2b, 3a, and 3b. Figures 2a and 2b show the 

15 images of Figures la and lb, respectively, after being edge-enhanced filtered. Figures 
3 a and 3b show the images of Figures 2a and 2b, respectively, after thresholding. The 
images of Figures 2a and 3 a show red channel intensity while the images of Figures 2b 
and 3b show infi-ared channel intensity. The images of Figures 2a and 2b show the 
result of filtering with an edge-enhancing filter to sharpen the boundaries between 

20 regions. The images of Figures 3a and 3b show binary images obtained by 

thresholding the images of Figures 2a and 2b. The white pixels in the images of 
Figures 3a and 3b indicate the presence of a sharp edge or high intensity •in the original 
image at the corresponding pixel. 

After edge-enhancement by filtering, the images are thresholded to 
25 create a binary image using a threshold for which N percent of the pixels are over the 
threshold. A value N=80 is preferred since it produced the best overall results. All 
pixels are thresholded to black or white independently, and there is no attempt to 
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create continuous lines. The threshold level in this image was chosen to cause 80% of 
the image pixels to survive thresholding and appear as white pixels. 

The prefiltering processing higUights both edges and low intensity 
regions in the images of Figures la and lb. The edges between light and dark and 

5 regions in the images of Figures 3 a and 3b generally correspond to edges in the 
original images of Figures la and lb, but may be offset slightly due to the action of 
the filter. The dark regions in the images of Figures 3a and 3b tend to be irregular and 
broken. There are two important characteristics of the image pair in Figures 3a and 
3b. Firstly, not all black regions are common to both images, and secondly, some of 

1 0 the ones in common have different fine structure. These observations indicate that 
correlations of the processed images will tend to have normalized correlation peaks 
below the ideal value of unity. These characteristics are not due to the specific edge- 
enhancement and filtering used, but are in fact due to the underlying differences in the 
original images of Figures la and lb. 

1 5 Optionally, the images of Figures 3a and 3b are blurred with a filter 

that thickens each pixel in the images by any means known in the art. One such way 
is to blur by t pixels in a square array that extends t pixels in all four directions from 
the central pixel Blurring may not be needed at full precision, but is preferable at 
reduced precision. 

20 After thresholding and blurring the images of Figures 3a and 3b, the 

resolution of the images are reduced by any means known in the art. One such way to 
reduce resolution by a factor of 2^^ is to partition the image into square blocks of pixels 
with 2^ pixels per side and to replace each square with the sum of the pixel values. 
This is equivalent to computing the low-low subband of the Haar wavelet of the image 

25 at level /. 
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The processed images are then registered. The resolution is then 
refined to obtain higher precision. Preferably, normalized correlation is used as the 
criterion for registration. 

The filtering, thresholding, blurring, resolution reduction, and 
5 registration described briefly above, will now be described in detail with reference to 
the Figures. The purpose of these steps in the registration methods of the present 
invention is to overcome the obstacles to successful registration mentioned earlier 
while retaining computational efficiency. 

The methods of the present invention use edge-enhanced images in 
1 0 order to capture information in edges, which is more rehable than pure intensity for 
multispectral images. However, it has been found experimentally that it is 
significantly better to threshold edge-enhanced images than to threshold edge-only 
images. Edge-enhancement creates strong regional boundaries, which tends to 
produce broader boundaries after thresholding than does edge-detection. Broad 
1 5 boundaries correlate better than do narrow boundaries. 

Many edge detectors use first-order directional derivatives to find 
transitions between regions. These detectors require multiple filter passes, each in a 
different primary direction. The registration methods of the present invention reduces 
the computational costs by eliminating the directional dependence, preferably by using 
20 second-derivative, direction-independent filters derived fi-om 2D Laplacian filters. 
The typical filter H has the form 



H 



-1 -1 "1 

-1 c -1 

-1 -1 -1 



(1) 



where C is a variable parameter. Because the result of filtering can be 
negative, the registration methods of the present invention use the absolute value of 
25 the filtered value rather than the signed value. A value C=8 creates an edge-only 
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filter, and sharp edges in the original appear as a pair of peaks in the filtered image. 
Values of C greater than 8 combine the edges with the image itself in different 
proportions, and thereby create an edge enhancement. 

Consider again the wetland scene depicted in the images of Figures la 
5 and lb. Notice how the natural features in the image of Figure 1 a differs from the 
natural features in the image of Figure lb. The infrared image of Figure lb reveals 
regions in the tideland that have almost uniform intensity in the infrared spectrum, but 
vary considerably in the visible spectrum. Manmade structures in the image of Figure 
lb tend to have similar boundaries in the two images. Both correlation and mutual 
1 0 information criteria tend to work well with the manmade structures but do poorly with 
the tidelands. Mutual information fails to register these images correctly, but 
normalized correlation of the unfiltered images succeeds, albeit with a low correlation 
value of 0.27. Both methods do poorly, in general, for this type of image, and in many 
cases, both methods fail. 

1 5 The edge-enhanced versions in Figures 2a and 2b corresponding to the 

images of Figures la and lb use the 9-point filter with a center value of 8.5. Note how 
much sharper the images are in Figures 2a and 2b, than they are in Figures la and lb. 

After filtering, the registration methods of the present invention 
thresholds the images of Figures 2a and 2b, preferably, to binary values. Because the 

20 images of Figures 2a and 2b contain some intensity information of the original images 
of Figures la and lb, the edges around regions of low intensity are less likely to 
survive thresholding than are the edges around high intensity regions. The threshold 
value is set in such a way as to pick up those edges, as well as some portions of higher 
intensity regions. The registration methods of the present invention preferably select 

25 the threshold automatically by computing a histogram of intensities and setting a 

threshold for which N% are over the threshold for a predetermined value of N. It has 
been found experimentally that the most effective thresholds for the filters used were 
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those for which 70 to 80% of the pixels were greater than the threshold. The images 
of Figures 3a and 3b use a threshold of 80%. 

The determination of the threshold is very similar to the histogram 
computation required to compute mutual information. The number of operations 

5 required is linearly proportional to the size of the image. The threshold computation is 
done just once per image registration, whereas the mutual information methods require 
the computation to be done for each relative position of the two images examined by 
the method. It is this aspect of the mutual information methods that incurs a heavy 
computational cost, and forces practical registration methods to do as few pairwise 

1 0 comparisons as possible. 

As an example of image characteristics that hinder registration, note the 
major regions of the images of Figures 3a and 3b that exhibit intensity reversal. Note 
also that the edges of the manmade structures are dark in both of the images of Figures 
3a and 3b. The nonuniformity of intensity reversal across the image greatly reduces 
1 5 the correlation peak heights, and decreases the likeUhood of a successful registration. 

Thickening (sometimes referred to as blurring) has two effects on 
correlation. The first is to enlarge the smaller features of an image so that they survive 
the filter/downsample process. This tends to increase the height of the correlation 
peak. The second effect is to broaden the correlation peak, which reduces the 
20 precision of the registration. Hence, thickening is useful to ensure that one can locate 
the correct registration position in a low-resolution image. Thickening is preferably 
abandoned or diminished at higher resolution in order to increase the precision of the 
final registration. 

Thickening is used to reduce the translation sensitivity of wavelet 
25 coefficients. The coefficients of a wavelet representation depend on the relative 
position of the image with respect to the underlying wavelet grid. A mathematical 
model and detailed experiments known in the art show that correlations of wavelet 
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coefficients in the low/low subband are relatively insensitive to translations, even 
though the wavelet coefficients themselves may be sensitive. However, this holds 
only for features large enough to be visible in the low-resolution wavelet subband. 
Note that for a resolution reduction of blocks of size 2' by 2' map into a single 

5 wavelet coefficient. Hence features of size on the order of 2' by 2' or smaller are too 
small to be captured well in the low/low subband coefficients, and they have little 
influence on the correlations of the wavelet coefficients. The idea behind thickening 
is to transform small features into larger ones that will be visible in the low/low 
wavelet subband. Enlarging these features tends to increase their participation in the 

1 0 correlation process. 

The choice of the thickening factor of t depends on the resolution of the 
wavelet. For example, for a resolution reduction of 1/64, blocks are size 8 by 8. 
Choosing t=\l ensures that features as small as 1 pixel will be visible in the low- 
resolution subband of the thickened image. However, correlation peaks broaden as 

1 5 thickening increases, making it more difficult to find the precise position of the 

correlation peak. One obtains better results overall by choosing a smaller value of t, 
say 5 or 9, for a resolution of 1/64. This compromise misses the smallest features in 
the images, which lowers the potential height of the correlation peak, but has littie 
impact on the correlation peak width. Thickening by t involves forming the sum of 

20 pixels of overlapping blocks of size 2^+ 1 x 2r+ 1 . It can be done efficientiy by a block- 
update calculation that scans the image fi-om left to right and top to bottom. The 
update requires only four operations per pixel plus a small overhead that depends on 
block size, but does not depend on the size of the image. 

Figures 4a, 4b, 5a, 5b, 6a, and 6b show thickening of 5, 9, and 17, 
25 respectively, of the images of Figures 3a and 3b. Note how thickening fills in the 

boundary Unes, Lines that are broken dots in the images of Figures 3a and 3b tend to 
be blurred solid lines in the corresponding images of Figures 4a, 4b, 5a, 5b, 6a, and 
6b. All of the images in Figures 4a, 4b, 5a, 5b, 6a, and 6b are shown at full resolution. 
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Figures 7a and 7b show the effect of resolution reduction after 
thickening of the image of Figure 3 a. The image of Figure 7a is the l/64th resolution 
reduction of the image of Figure 3 a, and the image of Figure 7b is the same for the 
image of Figure 6a (i.e., after thickening). Note that the vertical lines in the upper 
5 right of the image of Figure 7a are broken and imperfect because of the translation 
dependence of the downsampling and filtering in computing the Haar wavelet. In the 
image of Figure 7b, the same lines are more uniform because the image was blurred 
prior to computing wavelet coefficients. 

The filtering, thresholding, thickening, and wavelet subband operations 
10 can be done very efficiently. The process requires 1 0 floating-point operations per 
pixel to evaluate H, two to compute the histogram for the threshold, one for 
threshholding, four for thickening, and one for the Haar wavelet subband. This is a 
total of less than 20 floating-point operations per full-resolution pixel. Normalizations 
of these operations are not required because the normaUzed correlation coefficient 
15 calculation does all the normalization that is necessary. These filtering operations are 
done only once per registration. 

The computationally intensive part of a registration process is the 
evaluation of the registration criterion as a fimction of relative image position. If the 
cost is high, one must seek ways to keep the cost low. An effective way to do this 

20 involves a combination of resolution reduction and iteration. The registration methods 
of the present invention preferably use Fourier techniques in place of iteration to 
search large regions of the registration space very efficiently. This approach can be 
used effectively with resolution reduction to maintain low computational complexity. 
The key idea is that the normahzed correlation coefficient as a fimction of relative 

25 translational position reduces to a function of vector correlations. A brief summary 
for one-dimensional functions is the following. Let x = (xq, xi, . . . , xn-i) by an N- 
vector image, and y = (yo, yi, . . yM-i) be an M-vector pattern drawn fi:om a second 
image, with M < N. Let C(x, y)/ be the normalized correlation coefficient of y aligned 
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with X by shifting x relative to y by / positions, 0 < / <N-M. The normalized 
correlation in summation form is given by: 



C{x y) = - ^ r • (2) 

Equation (2) can be evaluated for all translations i for a cost equal to a 
5 small factor times the cost to evaluate it for a single value of /. The trick is to rewrite 

the equation in terms of vector correlations u O v defined to be 



(uG)v)i= Z«,^.v, (3) 

k 



where u and v are N vectors, and index expressions i+j are modulo N. 
Thus, four N- vectors are needed to convert Equation (2) into an equation that involves 

10 vector correlations in place of the summations that depend on /. Specifically, x is 
needed, and y is needed to extend to length N by appending N-M Os. This vector is 
denoted as ? . The vector whose elements are squares of the elements of x is also 
needed, which is denoted as x*^^ Finally, a mask vector m is needed whose first M 
elements are 1, and whose last N-M elements are 0. The mask m indicates which 

1 5 elements of y participate in the sums in Equation (2). In vector correlation form, 
Equation (2) becomes 

C{x,y)i= , ^ . ^ (4) 

/((x(2) O m), - (x O m)f ) [Et-o Vl " H {^'=0' V^)' 



Note that the summations of yk and yk^ in Equation (4) are independent 
of i and can be evaluated once per registration instead of once per relative position of 
20 the images. All N components of the vector correlations can be computed in the 
Fourier domain in a time proportional to N log N using fast Fourier transforms. 
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Equation (4) requires four Fourier transforms of real vectors and three inverse Fourier 
transforms to real vectors. Noting that a pair of real transforms can be performed 
forward and inversely as a single complex transform, the total cost is equal to two 
forward Fourier transforms and two inverse Fourier transforms of complex data. 

5 If N is not a power of 2, which is preferred for Fourier transforms, x 

can be extended to the next highest power of 2, with a corresponding mask vector for 
X. This results in a modified form of Equation (4) in which the summations of yk and 
yk^ become vector correlations involving the x mask, y and y . 

A Fourier-based registration search of all possible relative translations 
10 at any resolution can be done with about 500 to 700 flops (floating-point operations) 
per pixel at that resolution. This does not count the other operations per pixel or the 
fixed overhead in setting up the computation. Note that at 1/1 6th resolution, this is 
equivalent to about 33 to 45 flops per full resolution pixel, which is about twice the 
preprocessing cost. A mutual-information based registration method must do roughly 
15 20 to 40 operations per pixel for each relative position of images. Hence, the Fourier- 
based algorithm can examine all relative translations at a given precision for a cost 
equal to that incurred by mutual-information algorithms to examine a few dozen 
relative translations at the same precision. 

Iterative techniques usually require the estimate of a Hessian, gradient, 
20 and the mutual information function (or other criterion) to guide the direction of the 
search. Hence, the criterion function and its first and second derivatives must be well 
behaved over the search region in order to give correct estimates of the direction to 
move to improve registration. Because the Fourier-based method evaluates the 
registration criterion everywhere, it does not need to evaluate first and second 
25 derivatives, and is immune to problems caused by their potential ill behavior. 

The Fourier-based search described herein provides an alternative to 
iterative searches, and it may be much faster for some combinations of parameters. 
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Ultimately, its utility depends on the quality of its registrations. It is shown below that 
the Fourier-based correlation registration methods of the present invention produce 
results competitive with pixel-based mutual-information registration. 

EXAMPLES 

5 The experimental data reported below are the results of registrations of 

340 image sets whose characteristics are now described. Each data set contains four 
unregistered images ftom different bands. The three bands from the visible spectrum 
are red, blue, and green, and the fourth is infrared. 

These data sets were derived from 68 4-channel aerial image sets, each 
10 image of size 1536 by 1024. From each of these sets, five sets of size 512 by 512 
were extracted. Four sets were taken from the comers of the image and one set from 
the center. The misregistration from channel to channel was ±5 pixels in translation. 
Images differed as well by a very small rotation (a fraction of a degree) and by a small 
scale change. Within a 512 by 512 subimage, the rotation and scale change had little 
1 5 effect on the correlation. However, at the scale of the 1 536 by 1 024 image, the scale 
change and rotation were detectable and measurable. The translation offsets, scale 
change, and rotation values of the frill image from the registrations of the five 
subimages were able to be computed. 

The 68 frill images were drawn from four classes — agriculture (15 
20 images), forestry { 1 0 images), urban (33 images), and wetlands (1 0 images). Sample 
images of the first three appear in Figures 8a, 8b, 9a, 9b, 10a, and 10b, respectively, 
and the wetlands sample appears in Figures la and lb. As a class, the urban images 
tend to be the easiest to register because of the presence of sharp edges and comers 
that are visible across the spectra. In increasing order of difficulty are forestry, 
25 agriculture, and wetlands. Regisfration errors in the agricultural images are largely 
due to misregistration of similar features. For example, sfraight lines without 
crossings are very difficult to register. The agricultural and wetlands images tended to 
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be much more difficult because they contained fewer features in common across the 
spectra. 

There is no ground truth available for these data sets. However, sets of 
images within a class have channel offsets that are approximately equal for all images 
5 in that class, A consistent ground truth for all images of a class for one type of image 
were able to be computed by using this information. Some subimages in the wetlands 
set are essentially featureless because they are totally filled by a mudflat or water, and 
are impossible to register. 

The results for correlation-based registrations will now be presented 
10 without prefiltering, for pixel-based mutual information, and for correlation-based 
with prefiltering. In practice, the result of an image registration operation is a 
coordinate pair together with a number that represents the quality of that registration. 
For normalized correlation, the number is the height of the correlation peak. For 
mutual information, the number is the maximum of the mutual information function. 
15 If the quality measure is lower than a decision threshold, the registration is rejected. If 
the measure is equal or higher, the registration is considered valid, and the registration 
position is the position of the peak in the criterion function. 

When doing a registration for the mutual information criterion, a region 
of size 1 1 by 1 1 centered at ground truth was searched. Mutual information is 

20 prohibitively expensive if you search a large region exhaustively, and was very 
expensive even for the relatively small region that was searched. A more efficient 
approach would be to do an iterative search for the function maximum, however, the 
code is more complex and could be sensitive to the shape of the mutual information 
function. Instead, a complete search of a small local region centered at the correct 

25 registration was chosen. 

Normalized correlations were measured by using Fourier methods to 
build the global normalized correlation function as a function of relative position. 
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Within this function, a local search was conducted over an 1 1 by 11 square region 
centered at the correct registration, the same region over which a local search was 
conducted for the mutual-information-based search. Since 500 to 700 floating-point 
operations are required to search all 121 positions in the search space, this works out 
5 to 4 to 6 floating-point operations per pixel per position examined. This accounts for 
the low amortized cost of the Fourier-based search. 

The edge-enhancement step was tested on 5 different filters, four levels 
of thickening (1,5,9,17), and 5 levels of threshold (50%, 60%, 70%, 80%. and 90% of 
pixels over the threshold) for each of three resolutions. The filters used were 9-point 

10 filters with C=8, 8.5, 9.5, 10.5, and a 5-point filter with C=4. 100 parameter sets were 
applied at three resolutions to 1020 image pairs — 340 each of red-to-blue, red-to- 
green, and red-to-infrared registrations. This produced a total of 306,000 image 
registrations. In addition, pixel-based mutual-information registrations were 
performed at three resolutions, and correlation-based registrations of raw images at 

1 5 tiiree resolutions. Space restiictions limit this summary to the important highlights. 

Two measures — Recall and Precision — were used to evaluate the 
registrations. Each test produces a registration if the quality measure is over tiie 
decision threshold, and otherwise produces no registi-ation. A registiation was deemed 
correct if it matched ground truth +2 pixels. Hence there are three possible outcomes - 

20 ~ no match, correct registration, and false registration. The first measure. Recall, is 
the percentage of correct registrations out of the total number of images in the class. 
The second measure. Precision, is the percentage of correct registrations out of the 
sum of correct and incorrect registrations. Figure 1 1 contains plots known as receiver- 
operating characteristic (ROC) curves, which show the relation between Precision and 

25 Recall. Each point is a Recall/Precision pair for a particular setting of the registration- 
decision threshold. 

The ROC curves in Figure 1 1 compare the relative performance of 
registi-ation using correlation of filtered and unfiltered images with mutual- 



-19- 



information-based registration. The first column on the left shows the ROC curves for 
all images and all cross-spectral cases. Resolution goes from fcill resolution at the top 
row of Figure II, % resolution at the middle row, and to l/16th resolution at the 
bottom row. These curves illustrate that correlation of unfiUered images gives 
5 substantially poorer behavior over all the cases than do either normalized correlation 
of filtered images or maximization of mumal information. The filter chosen for this 
study is described in more detail below. Figure 1 1 reveals that filtering is effective in 
bringing normalized correlation to the point where it produces registrations 
comparable in quality to mutual information registrations. However, within the data 
1 0 set are subsets of images that are relative easy to register and some that are relatively 
difficult to register. The curves for the full data depend on the mixtures of those 
subsets in the full data set, and successes tend to mask failures. 

The second and third columns from the left of Figure 1 1 break up the 
data set into two subsets so that the performance on those subsets can be viewed 

1 5 individually. The second column plots ROC curves for red-to-infrared regisfrations, 
and they are clearly poorer than the comparable curves for the fiiU dataset in the first 
column. The third column plots the registrations of red-to-blue and red-to-green data, 
and they are clearly better than those in the first column. In fact, correlations of 
unfiltered images do very well on these images. This indicates that cross-spectral 

20 normalized correlation works well without filtering for band-to-band correlations in 
the visible spectrum. 

Figure 1 1 shows that the main difficulty for this data set Ues in the red- 
to-infrared registrations in the second column. Correlations of unfiltered images give 
very poor resuhs, but fihering brings the quality up to a range comparable to the 
25 quality of mutual information. 

In general, it is expected that Precision will fall with increasing Recall, 
because as the decision threshold is lowered to accept more registrations, a higher 
false-registration rate is likely. The curve in the upper left comer of Figure 1 1 has low 
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Precision at low Recall, and Precision increases with increasing Recall, which is 
unusual This occurs when, at a high registration-decision threshold, there are very 
few registrations accepted, most of which are incorrect. Hence both Precision and 
Recall are low. As the decision threshold decreases, more registrations are accepted, 
5 which boosts Recall, and if most of the registration decisions are correct. Precision 
also increases. 

Figures 12a and 12b show another way to partition the data set to 
illustrate the behavior of the registration methods. This partition is by type of image. 
The four graphs in Figure 12a are fiill resolution comparisons of urban, agriculture, 

10 forestry, and wetlands images for red to blue and green. All of the ROC curves for 
registration of red to blue and green are satisfactory for unfiltered correlation and 
mutual information registrations. Filtered correlation has some problems with the 
wetlands data, mainly because of featureless images. For cross-spectral registrations 
between visible spectra, Figure 12a indicates that normalized correlation of unfiltered 

1 5 images is satisfactory. 

Figure 12b contains similar data for registration of red to infrared, and 
reveals some difficult cases. It shows that unfiltered correlation performs poorly on 
the nonurban classes, and is not a viable approach for those images. Note that both 
maximization of mutual information and normalized correlation of filtered images 

20 perform well on these classes, with mutual information doing better on the agricultural 
images, and correlation of filtered images doing better on forestry and wetlands 
images. Agricultural images seem to be a problem for correlation of filtered images 
because regular field patterns often have multiple correlation peaks. Infrared 
intensities apparently lead some false peaks to be emphasized over the correct peaks. 

25 Lack of space does not permit us to show that performance falls off as resolution 
diminishes, and is otherwise consistent with the fixU resolution data. 

The main challenges for registering red to infrared are the agricultural 
and wetiands images. The search for suitable filters led to the choice of a filter with 
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edge-enhancement parameter C =8.5, edge threshold set to 80%, and no thickening. 
Figures 13 through 15 show the effect of varying filter parameters around this test set. 
In general, the parameter settings are robust in the sense that small changes of the 
parameters have only a small effect on performance. Also, no filter gives the best 
5 performance on all image classes and all cross-spectral cases. 

Figure 1 3 shows how performance varies with the choice of filter 
constant C. The two columns show the ROC curves for agricultural and wetlands 
images, respectively, and the three rows show Ml, 1/4, and l/16th resolution. All 
registrations are of red to infirared. For these plots, the filter for C=8, a pure 8-point 
1 0 Laplacian, has an edge-detection threshold set to 50%, for which it has its best 
performance. 

The effect of edge thresholding is illustrated in Figure 14. This figure 
is similar to Figure 13, except that the edge-detection threshold varies fi-om 50% to 
90% in each subplot. In all cases, the filter center value is C=8.5, and there is no edge 
1 5 thickening after edge detection. In this set of plots, the tiiresholds of 70% and 80% 
give similar performance. 

Figure 15 shows the effects of thickening on Recall and Precision. For 
wetlands data, thickening did not improve the registiration process, although its 
registration performance was very close to being equal to the unthickened data. 
20 Thickening with ? = 5 was slightly better than the unthickened data for the agriculhiral 
images. The data does not show that thickening helps at low resolution, as was 
expected would be the case. 

Note that the maximum recall rate for wetiands and agricultural data is 
on the order of 0.90 at fiiU resolution, and drops as resolution falls. For tiie wetiands 
25 data, about 1 0% of the images are unregisterable by ahnost any blind method for lack 
of common features. Figure 12 shows that, at full precision, mutual-information- 
based registi-ation was able to register most of agricultural images tiiat were not 
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registerable by normalized correlation. However, it had a lower maximum recall rate 
for the wetlands images. 

The combination of edge-enhancement, edge-detection, and Fourier- 
based normalized correlation is able to register images about as well as mutual- 
5 information-based methods, and is potentially faster. Fourier methods eliminate the 
need to use a nonlinear iteration to search for the relative translation that produces the 
best registration. The specific preprocessing steps investigated here appear to work 
well for cross-spectral registration of infrared to visible spectra, and may work across 
other spectra, provided that the images share a sufficient number of common features. 

10 Those skilled in the art will appreciate that the experiments discussed 

above, performed on a moderate-sized database, show that the registration methods of 
the present invention produced a correct registration rate of over 90% at a false 
positive rate of less than 10%. For a particularly difficult subset of images in the 
database, the correct registration rate fell to approximately 85% at a false positive rate 

1 5 of less than 1 0%. This retrieval quality is comparable to that of the mutual- 
information-based registration methods of the prior art. 

The methods of the present invention are particularly suited to be 
carried out by a computer software program, such computer software program 
preferably containing modules corresponding to the individual steps of the methods. 
20 Such software can of course be embodied in a computer-readable medium, such as an 
integrated chip or a peripheral device. 

While there has been shown and described what is considered to be 
preferred embodiments of the invention, it will, of course, be understood that various 
modifications and changes in form or detail could readily be made without departing 
25 from the spirit of the invention. It is therefore intended that the invention be not 

limited to the exact forms described and illustrated, but should be constructed to cover 
all modifications that may fall within the scope of the appended claims. 
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